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Abstract. Quantum Monte Carlo simulations of a two-component Bose mixture of 
trapped dipolar atoms of identical masses and dipole moments, provide numerical 
evidence of de-mixing at low finite temperatures. De-mixing occurs as a consequence 
of quantum statistics, which results in an effective attraction between like bosons. 
Spatial separation of two components takes place at low temperature with the onset 
of long exchanges of identical particles, underlying Bosc-Einstcin condensation of both 
components. Conversely, at higher temperature the system is miscible due to the 
entropy of mixing. Exchanges are also found to enhance demixing in the case of 
mixtures of non-identical and distinguishable species. 

PACS numbers: 67.60.Bc,67.85.-d,67.85.Hj,03.75.Mn 
1. Introduction 

Phase separation or de-mixing in multi-component mixtures has been a long standing 
topic of interest in chemistry and physics. Closely following the achievement of Bose- 
Einstein condensation in dilute gases [U [2] , there has been considerable interest in the 
study of binary mixtures of Bose-Einstein condensates (BECs). Of particular note, a 
two-component BEC was first reported in 1997 comprised of two hyperfine states of Rb 
[3], and then in 2001 using different atomic species (K and Rb) jlj. The advantage of 
these ultracold systems is that the entropy of mixing is small and de-mixing may be 
easily observed. 

The conditions under which de-mixing occurs in binary BEC mixtures with hard- 
core repulsion, have been the focus of a number of theoretical works, including mean- field 
treatments at zero [51 El [7J El [91 [10] and finite temperature [TTJ as well as Quantum Monte 
Carlo (QMC) simulations (121 [13]. Separation of species 1 and 2 is usually characterized 
in terms of a parameter A = UwU^i ~ U± 2 , defined in terms of the relative intraspecies 
(Un,U22) and interspecies (U12) interaction strengths. When A < 0, so that particles 
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of species 1 and 2 have a relatively strong repulsion, the system is predicted to phase 
separate, whereas for A > the system should remain mixed [HI El [7J El [91 [101 [TTJ [121 EE] • 
Recently, this criterion has been verified in experiments with binary BECs [H]. Earlier 
work on bosonic mixtures in the context of superfluid Helium have also predicted phase 
separation in the zero temperature limit for isotopes of different masses or concentrations 
[151 H6] - Indeed, all predictions of phase separation in binary mixtures of bosons to date 
rely on a mismatch of interaction strengths or of some other physical parameters on 
which the Hamiltonian depends (such as different particle masses, concentrations, or 
external trapping potentials for each species). We refer to this scenario as de-mixing 
through interactions. 

In this article, we report the prediction of de- mixing in a binary mixture of bosons 
with identical masses and interactions, due a very different mechanism — namely the 
effective attraction between indistinguishable bosons originating purely from quantum 
statistics. This prediction is made for A = in contrast to earlier work for which this 
value would lead to a miscible system. We refer to this scenario as de-mixing through 
exchanges or quantum de-mixing, which occurs when the system kinetic energy is re- 
duced by the formation of long exchanges of identical particles, leading to the spatial 
separation of the two components. This is, of course, the same mechanism underlying 
Bose-Einstein condensation, as first established by Feynman [T7J and subsequently elab- 
orated on [IH1 [19]. At low temperature, as the thermal wavelength becomes comparable 
to the interparticle distance, quantum exchanges involving two or more indistinguishable 
particles become frequent, and condensation sets in - this effect is also responsible for 
phase separation. 

For this study, we have elected to use the dipole-dipole interaction potential to 
describe the inter- and intra-species interactions, for which the condition A = is al- 
ways satisfied. The physics of ultracold dipolar bosons has fast become the subject of 
intense research activity. Bose-Einstein condensation of dipolar Chromium atoms has 
already been achieved [20]. The long-range and anisotropic nature of the interaction 
leads to many fascinating and novel phenomena (see the review [21] and references con- 
tained therein). 

So far, there have been few calculations explicitly dealing with binary mixtures 
of dipolar bosons j22j [23]. In Ref. [22], the stability of a binary mixture with the 
components having oppositely oriented dipoles was investigated. In Ref. [23], sponta- 
neous pattern formation associated with ferrofluidity was predicted and attributed to 
the anisotropic nature of the interaction. The role of finite range interactions in binary 
mixtures has also been addressed in previous studies [101 E3] - In particular, it was found 
that increasing the range of the interactions leads to increased mixing |10j . 
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2. Formulation 



We consider here a system comprising N a atoms of species a and iV& atoms of species b 
confined in a harmonic trapping potential. Let M a , Mb be the masses of each species, 
and V aa , Vbb and V a b be the inter- and intra-species interaction potentials. We consider 
here a two-dimensional confined geometry, with a transverse polarization field in the z 
direction,for which the dipole-dipole interaction potential becomes isotropic and purely 
repulsive up, i.e., V mrn i(r) = d m d' m /r 3 , between particles of species m and m', at a distance 
r from each other and with respective (electric or magnetic) dipole moments d m and d' m . 
Introducing a reference dipole moment d , we choose characteristic length and energy 
scales as r D = d 2 Q M a /h 2 and e Q = h 2 /(M a r 2 ) respectively. The Hamiltonian for the 
two-component system is given, in dimensionless form, by: 




where r mk is the position of the kth particle of species m and T = l/2(L p /r )~ 4 gives 
the trap strength, L p = ^Jh/M a uj being the harmonic oscillator length. For brevity, we 
write the relative dipole amplitudes as a = d a /d and = db/d . 

If by analogy with the case of hardcore boson mixtures, we define a mixing 
parameter A ~ V aa V aa — V 2 b , then A = always. 



3. Numerical results 



Henceforth, we choose M a = Mb and N a = Nb = N/2. We have investigated the finite 
temperature equilibrium properties of the system by QMC simulations based on the 
Continuous-space Worm Algorithm [25| [26] . This technique is numerically exact, to 
within a controllable statistical error. We have carried out simulations with two values 
of N, namely 40 and 100, and with different values of the harmonic trap strength T, 
always chosen sufficiently small to keep the density in the middle of the trap below 
the crystallization threshold (27]. Our results are qualitatively the same for all cases 
considered. 



3.1. Identical and distinguishable species 

We initially set the dipole moments of each species equal (a = (3 = 1) so that the 
Hamiltonian (P) is symmetric to an interchange of labels a and b. Typical QMC 



| This is in fact a valid approximation for so-called pancake traps — highly anisotropic harmonic traps 
with transverse and planar trapping frequencies satisfying lu z lo p . 




Figure 1. (Color online) QMC configuration snapshots at temperature T = 0.5 for 
species a and b in harmonic trap with a = /3 = 1. Snapshots show particle world lines. 
In the upper plots N = 40 and T = 8 whereas in the lower plots N = 100 and T = 0.5. 



configuration snapshots are shown in Fig. CD for N = 40 with T = 8 (top) and N = 100 
with T = 0.5 (bottom). In each case, the left plot shows the case where exchanges 
between indistinguishable (like) particles are included, whereas the right plot shows 
the same system but without exchanges (ie. distinguishable particle statistics). For 
indistinguishable particle statistics like particles tend to aggregate, which is suggestive 
of de-mixing, whereas for distinguishable particle statistics the system remains mixed. 
Individual snapshots are not conclusive however, so this prediction is verified first in 
Fig. |2] for the case N = AO and T = 8 and then in Fig. [3] for a larger system size with 
N = 100 and T = 0.5. Both figures show the integrated [§| pair correlation function 
gab{r), which for the lowest temperature shown, T = 0.5, becomes suppressed at short 
distances with respect to g aa { r )- That is, the probability of finding unlike particles 
separated by r is less than that of finding like particles at the same distance. Differently 
phrased, a particle of a given species is preferentially surrounded by like particles. For 
both system sizes, at the highest temperature shown, there is no evidence of de-mixing, 
and the probability density of position for particles of either species only depends on 
the distance from the center of the trap. In this regime, entropy dominates. 

§ Clearly, since the system is not translationally invariant one ought more properly look at the two- 
point correlation function y mm '(r, r'), with r, r' measured with respect to the center of the trap. The 
function considered here is averaged over the whole trap; however, it still provides the quantitative 
information sought here. 
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Figure 3. (Color online) Integrated pair correlation function at temperatures T = 
5, 1,0.5 between like (g aa ) and unlike (g a b) particles for system in harmonic trap with 
r = 0.5, a 2 = (3 2 = 1 and N = 100. Errors arc < 2 x 10" 2 in all cases. 
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Figure 4. Relative frequency of permutation cycles of length L for species a with 
N = 100, at four different temperatures for the case where a 2 = (5 2 = 1. 

As the temperature is lowered, the two components separate. To verify that 
exchanges are responsible for the observed de-mixing Figs. [2] (bottom) and [3] (bottom) 
also show the result for T=0.5, but with exchanges turned off in the simulation (i.e., 
we regard particles as distinguishable). The resulting plots shows the system is fully 
mixed. The energetic mechanism leading to phase segregation, is that particles can lower 
their kinetic energy by exchanging with like particles, thereby enhancing their spatial 
derealization. We also expect de-mixing to occur for larger N (relevant to experiments) 
as the relative interfacial energy decreases. 

It should be emphasized that the de-mixing predicted here is a finite temperature 
effect. Finally, on lowering the temperature even further, the system remixes. This 
effect, also observed in hardcore boson mixtures [28], occurs when the temperature T 
is less than the level spacing hu = v2T of the confining harmonic potential. We have 
verified in our simulations that the remixing temperature decreases for smaller T. 

We further elaborate on the connection between the appearance of long permutation 
cycles and de-mixing referring to Fig. HI which shows the relative frequency of 
permutation cycles including different numbers L of particles, for a single species. Note 
that, as T — > 0, the frequency of occurrence of cycles of permutation (ie. exchanges) 
involving almost all particles of each species (i.e., L = 50, in this case) increases 
dramatically. These exchanges are central to Bose- Einstein condensation, which in turn 
is responsible for the effective attraction between like atoms and therefore the observed 
phase separation. Such quantum de-mixing is not particular to our choice of interaction. 
Indeed, we have verified that it also occurs for potentials of the form V(r) ~ 1/r 12 (and 
with A = 0), which emulate the usual hard-core repulsion appropriate for ultracold 
alkali atoms. In the more general case where A > the repulsive interaction energy 
would be offset by the decrease of kinetic energy due to exchanges, although this effect 
may not be large enough to lead to phase demixing. 

It is instructive to note that only a technique which explicitly treats exchanges of 
indistinguishable particles at finite temperature, can yield predictions of de-mixing such 
as those shown here, for components of identical masses and interactions. 
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Figure 5. (Color online) Density profiles (left) and pair correlation functions (right) 
at temperatures T = 10, f, 0.5, for species a and b in harmonic trap with T = 0.5, 
a 2 = 1,(3 2 = 1.21 and N = 100. Also repeated for comparison is the case T = 0.5 
without exchanges (i.e., distinguishable quantum particles). Statistical errors are 
< 2 x 10~ 2 for all densities and < 5 x 10~ 2 for all pair correlation functions shown. 
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Figure 6. (Color online) Radial mixing A as a function of temperature T for: (a) 
a 2 = 1,(3 2 = 1.21; (b) a 2 = 0.8, (3 2 = 1.25; and (c) a 2 = 0.5, (3 2 = 2. 



3.2. Non-identical and distinguishable species 

Next, we consider the situation where the dipole moments of each species are not equal 
(ie. a 7^ /3). In this case, the symmetry of the Hamiltonian with respect of exchange of 
species a and b is explicitly broken, and we can expect phase separation to be evident 
in the radial density profiles p m (r), m = a,b, computed with respect to the center of 
the trap. 

We consider for definiteness a < (3. In this case, species b forms a shell around 
species a due to the higher interaction energy of species b in the presence of the 
harmonic trap. That is, the larger interparticle repulsion of species b particles pushes 
that component to the outside of the trap. This can be seen in Fig. |5]for the case where 
a 2 = 1,(3 2 = 1.21. Specifically, as the temperature is lowered the partial overlap in 
density profiles between species a and b decreases, indicating the onset of de-mixing. 
It is worth noting that for binary mixtures of non-identical species the presence of de- 
mixing is due to the combined effects of interactions and quantum exchanges. To verify 
that exchanges do indeed enhance de-mixing, we show in Fig. [5] density profiles and pair 
correlation functions at T = 0.5, for the simulations both with and without exchanges. 

The functions p m {r) can be used to quantify the degree of phase separation through 
the normalized overlap integral: 



When there is complete overlap (p a oc pb) then A = 1 indicating total mixing, whereas 
for complete phase separation we have A = 0. Note that A is not a good indicator 
of de-mixing when a = (3 since in this case the symmetry of the Hamiltonian means 
that the density profiles are identical when sufficient statistics are accumulated in the 
simulations. Fig. [6] shows the radial mixing A as a function of temperature for three 
different combinations of dipole strengths (a and (3). As the temperature decreases 
and/or the ratio between a and (3 decreases phase separation becomes more pronounced. 



A 



[/ pa(r)p b (r)dr] 



(2) 



[/ Pa(r) 2 dr] [f p b (r) 2 dr] 
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4. Conclusions 

In conclusion, low-temperature de-mixing purely due to Bose statistics is predicted to 
occur in a binary mixture of dipolar atoms, even when masses and the dipole moments 
of each species are equal, i.e., inter- and intra-species interactions are identical. In the 
case where the dipole moments are unequal, interactions also contribute to de-mixing. 
As the temperature is raised, on the other hand, the system becomes miscible due to 
the entropy of mixing. The observation of the effect predicted here appears well within 
reach of current experimental efforts with cold dipolar systems. For our simulations the 
total peak density was pr 2 ~ 1, which for 52 Cr atoms gives p ~ 10 17 m~ 2 . This, while 
higher than those typical of magnetic traps, should be approachable using optical traps 
[27]. Moreover, optical traps allow for the simultaneous trapping of different hyperfme 
states [29], which in principle allows for two distinct but equal mass species. The inte- 
grated pair correlation function should be measurable by averaging over several "single 
shot" absorption images [30] . 
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